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■■ ■ Abstract 
00 ■ 

The recently discovered A (3872) has many possible interpretations. We study the production 
of A(3872) with PANDA at GSI for the antiproton-proton collision with two possible interpre- 
tations of X(3872). One is as a loosely-bound molecule of D-mesons, while another is a 2P 
charmonium state Xci (2P). Using effective couplings we are able to give numerical predictions 
for the production near the threshold and the production associated with ir°. The produced 
A (3872) can be identified with its decay J/ipir + ir~. We also study the possible background 
£ — , ' near the threshold production for A(3872) — > J/ip-K + TT~ . With the designed luminosity 1.5fb _1 

per year of PANDA we find that the event number of pp — ► J/ipir + Tr~ near the threshold is 
at the order of 10 6 ~ 10 8 , where the large uncertainty comes from the total decay width of 
' (**]'■ A (3872). Our study shows that at the threshold more than about 60% events come from the 

<*S decay of A(3872) and two interpretations are distinguishable from the line-shape of the produc- 

JL , tion. With our results we except that the PANDA experiments will shed light on the property 

ofA(3872). 

> 

OS ■ A" (3872) has been first discovered by Belle collaboration [1] in the decays B — ► If X (3872) 

K J / ifjTr + tt^ . Later, its existence has been confirmed by experiments of Babar[2j, CDF[3] and 
D0@]. The word average mass of X(3872) now is m x = (3871.2 ± 0.5)MeV and the total width 
is Tx < 2.3MeV at 90% C.L.jS]. The angular distribution analysis made by Belle [6j favors 
J PC = 1++. A similar analysis by CDF [7] collaboration allow J PC = 1++ and J PC = 2~+ as 
well. The dipion mass distribution in the decay into J/tfjTr + ir~ suggests that the tt + tt~ may come 
from a p resonance, this is supported by the CDF analysis [7j. 

Many interpretations of X(3872) exist. In [8] it is interpreted as a loosely-bound molecule of 
D D — c.c In P [10] it is suggested that X(3872) is the first excited state of the conventional 
charmonium Xci> i-e., Xci(2P)- Other possible interpretations are also possible, like the /S-wave 
threshold effect of D°D*° [IT], a cusp effect ^2]. a diquark anti-diquark bound state[l_3], a hybrid 
charmonium state|14j and a tetraquark state|15j. etc. The existence of these many interpretations 
reflects the fact that the structure of AT (3872) is still unclear. It is clear that further studies in 
experiment and theory are needed. 

In this work we study the production of X(3872) in pp collisions by taking X(3872) as a loosely- 
bound molecule of D D*° — c.c. or as the first excited state of the conventional charmonium Xci- 
Experimentally the production can be studied with PANDA detector for the anti-proton beam 
facility at GSI[16j, where the anti-proton is with the energy from 1 ~ 15GeV. In pp collisions 
AT (3872) can be produced near its threshold. We assume it will be identified through its decay into 
J/ipTr + Tr~ , then the same final state can also be produced through direct production, which will 
be a background in identification of AT(3872). We will make numerical predictions for the process 
pp — > J/ipiT + ir~ near the threshold of A"(3872), where the final state is produced through the decay 
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of X(3872) or through the direct production. We will also give numerical results for the production 
associated with a tt°. Theoretical study of the X(3872) production at quark-gluon level in the 
energy range we consider is very difficult. We will take the approach of effective Lagrangian in 
terms of hadrons. We first discuss couplings between relevant hadrons and then give our numerical 
results. 

If we assume the X(3872) is a pure 2P charmonium state %& (2P), then we can estimate it decay 
width of into pp as following. In the decay the charm quark pair will be annihilated into gluons first, 
then those gluons will be converted into the pp pair. The conversion will be the same for Xd in the 
ground and the first excited state. We take charm quarks as heavy quarks and use a nonrelativistic 
wave function to describe the charm quark pair in the charmonia. In the nonrelativistic limit, the 
annihilation rate of Xd into gluons will be proportional to the square of the first derivative of the 
radial wave-function R(r) . Therefore we have: 

r[X(3872)-pp] _ I#(0)IL(2P)_ (1) 



nxd^pp] mom 
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One can obtain the wave functions with some potential models. In [TJj the numerical results for 
four different potentials are given. Here we use the result with the Cornell potential |17|: 



I^(°)IL(2P) _ 0.186 



#(o)IL °- 131 



1.42. (2) 



From other three models the ratio is 0.97, 1.05 and 1.33, respectively. One can re-scale our pre- 
diction for the total cross-section with the ratio from the Cornell model to obtain the prediction 
with ratios from other three models. It should be noted that in [T7] the main quantum number is 
defined as n r + 1 + 1. Therefore the 2P state in [T7] is the P-wave ground state while the 3P state 
is the first excited P-wave state. Using the above results and experimental data we can determine 
the effective coupling constant g p px which is defined as 

£ P px = gppxVl^lhpX^ g pP x = 1-11 x 10~ 3 , (3) 

where is the effective field for X(3872). 

If X(3872) is a loosely-bound molecule of D°D*°, the decay width into pp has been estimated 
by pi] as: 

r[X(3872)^pp]= (A^_^_y /2 (35e v), Ex = M D0D - t0 -mx, (4) 

where A can be chosen as m n since low-energy scattering of charm mesons is dominated by pion 
exchange and Ex = 0.6 ±0.6 MeV is the bounding energy of the molecular state. We use A = m n , 
Ex = 0.6 MeV and have for the effective coupling: 

g ppx = 7.14 x 10" 4 . (5) 

Having fixed the coupling with pp we turn to the decay X — > J/ipir + ir~. As discussed at 
the beginning, it is likely that the 7r-pair comes from the p-resonance. We will take the decay as 
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X — > J/ipp — ► J/vfjir + ir . Then decay amplitude with effective couplings can be written as: 



A v [p -> 7T + 7T~] = ^G p7T7T (p + - p-) u , 

A^a{X -> J/ipp] = Gxippe^apq^ef , (6) 

where q is the four momentum of p, p+ and p- are the momentum of 7r + ,7r~, respectively, and 
ex, are the polarization four vector of the X(3872) and the J/ip. The coupling constant G pnn 
can be determined from the decay of p into tt + tt~ , which is 11.99 ± 0.06. The total decay width 
can be obtained as: 

F[X - J/Vvr+vr-] = |G x ^ p | 2 (226keV), (7) 

where we have used a cutoff for the invariant mass of the -zr-pair, which is taken as mi-^ > 400MeV 
as in the experiment of Belle pQ. 

If X(3872) is a loosely-bound state of the charm mesons, the coupling Gxipp can be expressed 
with the total width and binding energy [2D]: 



'Ex + T 2 x /(16E x )\ 



1/2 



where T x is the total width of the X(3872) and the lower bound on width to be T x > 2T[D*°] = 
136 ± 32keV[2D]. If we take E x = 0.6MeV and the upper bound to be 2.3 MeV, then we obtain 

Gxipp ~ 0.893 ~ 1.05. (9) 

For the case that X(3872) is a 2P charmonium state Xci(2P), the decay width is estimated as[9j: 

F[X -» J/r/nr + n~] = 40keV, (10) 

which gives the value of the effective coupling: 

G x ^ p « 0.42. (11) 

With the estimated coupling constants between relevant hadrons we are able to predict the 
cross-section for the process pp — > J /ipir + ir~ near the threshold of X(3872), where the final state 
can be produced from the decay of X(3872) or from direct production. The amplitude for the final 
state from the decay can be written as: 

pp _ pp pv 12 

M \pp^X^ J/Vvr+vr-] = Ap[pp -> X] £ _ ^ - ,^ * A va [X -> J/#] 

W aP A , + - 

7T 7T 



q z — m~ + im p L p 

•Ap\pp^X] = gppxvpsjp^Ups, (12) 

where P denotes the momentum of the pp-pair, and q is the momentum of the 7r-pair. The final 
state can also be produced directly from the pp- annihilation as shown in Fig.l., this should be 
taken as a background for the production of X(3872). 
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Figure 1: The production of J/i[m + ir through the pp-annihilation 

In Fig.l n denotes the internal neutron line. Again we use effective couplings to calculate the 
process. We use — v2<?pp7r75 for the npir effective vertex, — <7 pp7r 75 for the ppir effective vertex, 
Wppiplii for the J/tppp, and ig n nii>lp, for the J/ipnn, respectively. The amplitude from Fig.f can be 
expressed as 

M a = x^g^Vprts _ j } _ mn 75 {f _ _ mp W^, 
M b = i2g^g m v pslfl ^ _ * _ - 75 ^ _ -L _ — 

-Mc = ^ n ^ ps - 75 _ 1 _ ^ ^ _ ^ _ — 75^6^, (13) 

then the total amplitude is the sum: 

M[pp -» J/^it + it~} = M[pp^ X -> J/Vot + 7t~] + M a + M b + M c . (14) 

The effective coupling g ppn is g ppn = 13.5. By using isospin symmetry we have <7 pp ^ = g n nip and 
m n = m p . The effective coupling g ppi p can be determined from the decay J/tp — > pp. It should 
be noted that for the decay it is possible that another coupling, i.e., the Pauli's coupling can be 
appear[18]. We neglect this coupling and get g ppi p = 1.62 x 1CP 3 [18J. Although the coupling 
constants are estimated, but their relative phase is unknown. There are two possible cases: Case 
1: The product gpp^gp 1 ^ has the same sign as that of the product g vp xGx^pG p7T7T . Case 2: The 
two products have different sign. The expression pf the amplitude squared is too length because it 
is a 2 — > 3 body process and we do not try to produce an analytical expression for the total cross 
section. Instead giving the analytical expression we simply take the amplitude squared to perform 
the phase space integral numerically. In Fig. 2 and Fig. 3 we plot the total cross section as functions 
of the invariant mass s of the pp pair, where we take Tx = 2.3 MeV and the coupling constants 
estimated before. Fig. 2 is for X(3872) as a loosely bound state of .D-mesons, Fig. 3 is for X(3872) 
as Xd(2P). 

From these figures we clearly see that from the line shape of the cross-section two interpretations 
of X(3872) can be distinguished. The difference comes from the decay into J/ip7r + ir~ with the 
different assignment of X(3872). We also see that the background is an significant contribution 
for the production. For the assignment with the bound state of D-mesons we have the total cross- 
section at y/a = 3.872GeV by taking E x = 0.6 MeV, T x = 136keV ~ 2.3MeV: 

<j[pp -> X -> J/ijm + n~] = 3.57 ~ 443nb, (15) 
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Figure 2: The total cross-section as a function of y/s for X(3872) as a loosely bound state of D- 
mesons. The lower resonance curve is without the background, the upper curves are with the back 
ground. The solid one is for Case 2, while the dashed one is for Case 1. 

for the assignment with Xci{^P) we have with Tx = 136keV ~ 2.3MeV: 

a[pp -» X -> J/V'vr+vr"] = 2.19 ~ 238nb. (16) 

The main uncertainties in the above come from the unknown width Tx- In Fig. 4 we plot the 
dependence of the cross section with different assignments as a function of Tx- From Fig. 4 we see 
that there is a strong dependence of the total cross section on the total width and the cross-section 
with two interpretations are different. Hence the measurement of the cross section will give a clear 
evidence to indicate which interpretation is the correct one. 

The luminosity of PANDA can be up to 2 x 10 32 cm _2 s _1 [16J. Assuming 50% overall efficiency 
and 6 months/year data taking, the integrated luminosity is to be 1.5fb _1 per year. With the 
integrated luminosity and the cross-section obtained here, one can expect 10 6 ~ 10 8 events per 
year for the production of J/ipTr + ir~ near the threshold. With the large number of events one can 
study X(3872) in more detail. 

With the estimated couplings we can also study the production of X(3872) associated with ir°, 
i.e., the production away from the resonance region. There are two diagrams for the process given 
in Fig. 5. The calculation of the total cross section is straightforward. The analytical expression for 
the differential cross-section can be found in [18] . We only give our numerical result here. With 
the same parameters we obtain the total cross section as a function of yfs up to 5GeV given in 
Fig. 6. From Fig. 6 we find that the total cross section of pp — > Xi: is at order of ~ lOOpb. With 
the designed luminosity and by considering the branching ratio of decays of X it is likely that such 
a process can not be observed. 

To summarize: In this work we have studied the A(3872)-production at PANDA, where two 
possible interpretations of A(3872) have been assumed. We have found that there will be a large 
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Figure 3: The total cross-section as a function of \fs for X(3872) as Xci(2-P). The lower resonance 
curve is without the background, the upper curves are with the back ground. The solid one is for 
Case 2, while the dashed one is for Case 1. 

number of events for the process pp — > J/^tt + tt^ at the threshold where large fraction of events 
will be produced from the decay of X(3872). By measuring the cross-section and its s-dependence 
near the threshold one can distinguish the two interpretations. For other possible interpretations 
like a diquark anti-diquark bound state, a hybrid charmonium state and a tetraquark state, the 
coupling with pp is so far unknown. Once the coupling is estimated, the production rate can 
be obtained from our results here. If the coupling is not extremely small in comparison with 
those given in Eq.(3,5), one may still expect that X(3872) can be produced with a not small event 
number. Hence the study of the AT(3872)-production at PANDA will provide important information 
about the structure of X(3872). We have also studied the production associated with tt°. But the 
cross-section by considering the branching ratio of X(3872) may be too small to be measured. 
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